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Abstract 

A model kinetic equation is constructed for the transport of a massless 
Bose gas. This equation is applied to solution of the boundary value prob- 
lem for the transport of radiation in the half-space occupied by a dispersive 
medium that is in local thermal equilibrium with the radiation. It is shown 
that the difference in temperature between the dispersive medium and the 
incident radiation depends substantially on the character of the scattering 
properties of the particles of the medium. 
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02.30.Rz Integral equations, 05. 70. -a Thermodynamics. 

Introduction. Statement of the problem and basic equations 

In recent decades, transport theory has been substantially devel- 
oped for Boltzmann gas [1, 2], Boltzmann plasma [3, 4], and Fermi gas 
[3-6]. Progress in describing radiation transport in dispersive ("tur- 
bid") media should also be mentioned [7, 8]. Both approximate and 
exact analytic methods (for different collision - integral models) have 
been used in these investigations. 

In contrast to the case of Boltzmann and Fermi gases, no general 
statement of the transport problem is known for Bose gas, although 
some particular results have been obtained in this case (e.g., see [7, 8]). 
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At the same time, transport phenomena are frequently encountered 
in practice, e.g., heat transport by phonons in a sohd and radiative 
transport by photons in dispersive media. All of these processes have 
a number of common properties. In what follows, we confine ourselves 
to the case of massless Bose systems for which the chemical potential 
vanishes in the thermodynamic equilibrium. 

Consider the problem of constructing a transport equation describ- 
ing the evolution of Bose gas. The collision integrals used in trans- 
port theory take into account either the elastic scattering of photons 
(electromagnetic waves) or their absorption by a dispersive medium. 
However, if the dispersive medium in which radiative transport takes 
place is in local equilibrium with the photon gas, then it not only ab- 
sorbs but also re-radiates the photons. As an example, we mention 
the case of radiative transport in gaseous dust clouds in the cosmos, 
where precisely this type of radiative transport from central regions 
to the outside of the cloud results in cooling and compression of the 
cloud, which, by-turn, gives rise to the star-formation process [9] . 

As an example of the construction of a model transport equation, 
we consider the Bhatnagar — Gross — Krook (BGK) model for the 
Boltzmann equation. The Boltzmann equation for the distribution 
function f{r,\,t) of the velocity of gas molecules has the form 



where J[f] is the collision integral. 

In constructing the model, the left-hand side of Eq. (1) remains 
unchanged, whereas the collision integral is transformed according to 
the following scheme. It is assumed that the molecules come to the 
equilibrium distribution after a single collision, i.e., the molecule dic- 
tistribution function after the collision coincides with the equilibrium 
distribution function 



dt 



+ vV/ = J[f] 
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(2) 



where ^{v) is the effective colhsion frequency. 

The Boltzmann equation (1) is compatible with the conservation 
laws for the number of molecules, momentum, and energy. The model 
equation (2) must possess similar properties. This requirement leads 
to the uniquely determined parameters no, Tq and u (for more detail, 
see [10]). In the case iy{v) = const, the parameters no,To and u 
coincide with the average molecule concentration, temperature, and 
velocity, respectively. 

We consider the case of a massless Bose gas, meaning a photon 
gas. The case of a phonon gas differs only in some details, which are 
indicated in the course of the presentation. 

The state of a photon gas is described by the distribution func- 
tion fi{uj, r, n), where oj is the frequency, n is the direction of photon 
propagation (n = — k, k is the wave vector and c is the velocity of 



light), and the subscript i (i = 1, 2) indexes the polarization direction. 
For a phonon gas, there are three possible types of polarization and, 
therefore, the subscript i assumes the values i = 1,2,3 in this case. 
The transport equation for the distribution function r, n) has the 
form 



Here the functional J[f] describes the absorption and re-radiation 
process for phonons in the dispersive medium where they propagate. 
We note that in contrast to photons, the phonons interact with one 
another. In this case, J[f] means the collision integral for phonons. 
In Eq. (3) , the parameter c is the velocity of light for the photon gas 
and the velocity of sound in the case of the phonon gas (on the condi- 
tion that the propagation velocities of the longitudinal and transverse 



CO 



-^ + cnWf = J[f\. 



(3) 
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oscillation coincide). 

Let a photon gas be in local equilibrium with a dispersive medium. 

In this case, the distribution function of photons re-radiated by the 
medium coincides with the Planck function 

r ~^ 

fp = 



fyjj 



where T is a parameter defined below. Equation (3) takes the form 

df 

-Al + cnVfi = uiuj)ifp-fi), (4) 

where z^(cj) is the effective froquency of the phonon collisions. For pho- 
tons, this quantity is related to the effective absorption cross-section 
for the dispersing particles, namely. 



z/(a;) ~ ^ J 'n'{r)(7r{oj)dr, 



where n(r) is the radial distribution of the particles (on the condition 
that they have a spherical form) and (7r{oj) is the absorption cross- 
section for a particle of radius r. 

Note that in the general case, it is necessary to distinguish between 
photon scattering and re-radiation. For simplicity, we do not discuss 
this question in the present paper. 

In general, the number of photons and their momenta are not pre- 
served under re-radiation. At the same time, the photon gas energy 
must be invariant in local thermodynamic equilibrium between the 
photon gas and the dispersion system. This means that the integral 
on the right-hand side of Eq. (4) vanishes, 

f Y,^M^)Up - fi)dnduj = 0. (5) 

i=l 

Relationship (5) is the one whichn determines the parameter T. 
We consider the case where the Bose system is in a nearly equilib- 
rium state. Here the distribution function fi can be represented in the 
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form 



(6) 



where ipi is a linear correction to the Planck distribution fp \ 

In what follows, we assume that the radiation is not polarized and 
the distribution function does not depend on i. Therefore, we omit 
the subscript i in the distribution functions and i/ji. The closeness 
to equilibrium means that Tq and the parameter T, which has the 
meaning of the effective temperature, vary weakly in the gas volume, 
which is equivalent to the relationship 



l^i « 1, 



(7) 



where S = — — — , and Tq is the temperature of the Bose gas at some 
point in the volume. 

Under the condition (7), we can linearize the Planck function fp{6) 
with respect to the parameter 6, fp{5) = /p(0) + /p(0)J, where 

-1 -1 



or, in explicit form. 



exp 



KTo + 6To) 



1 



Here E{x) is the Einstein function. 

On substituting (6)- (8) into Eq. (4), we derive 



(8) 



dijj 
'dt 



huj huj 



(9) 



Equation (9) is the linearized form of the transport equation (4) . 
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It follows from (5) that the difference T—To is given by the formula 

-1 



T-n = 



kTr 







kkTqJ 



iy{uj)uj^iljdujdn. 



Equation (9) shows that it is convenient to introduce a new function 
ip instead of tp, 

hjj \ 

and, then, Eq. (9) can be rewritten as 



ip{uj, r, n) = ujE(^^^(f{uj, r, n). 



dt 



+ cnV(p = 



where 



/uUjj')uj'^ E {^^^^^ipUjj\ r, n) du'dn — i'(u)ip(u, r, n), (10) 
\kToy 

1 = j uiu'^E^^^du'dn. 

We consider the case of a stationary process and assume that the 
function z/(cj) can be approximated by a power function, z/(cj) = uqU)'^. 
Note that the case a = 2 corresponds to the absorption of electromag- 
netic radiation by particles whose size is small in comparison with 
the wavelength. For large-sized particles, a = 0. We introduce the 
dimensionless variables 



UJ 



and r* = — 



kTo c \ h 

In this case, Eq. (10) takes the form (here and henceforth, the 
asterisks in the variables are omitted) 



u'^^^E{uj')ip{uj\ r, n) duj'dn — u'^ip^tu, r, n), (11) 



47r/o(c>^) 



where 



00 



47Tlo{a) = J uj''+^E{uj)dujdn = J ^'^^^E{u)du, 
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hence 

00 





We note that it is assumed in the above calculations that does 
not depend on the coordinate, though, generally, this is not the case. 
This fact can be taken into account by introducing the optical length 
[8]. We do not dwell on this in what follows. 

We consider the one-dimensional problem in which ip = (fi{x,v,ij), 
where v = n^jV, i.e., v is the cosine of the angle between the direction of 
the quantum motion and the x-axis. Here, Eq. (11) can be rewritten 
as 

d 

v—(f(x, V, cj) + cj^(f(x, V, cu) = 
ox 

00 1 



^ / /a+4 



2k{a) J 





E{uj')duj' j (p{x,v',ij')dv', (12) 



where 



00 00 





In the transition to the Eq. (12), we integrated with respect to the 
azimuthal angle from to 27r because the function cp does not depend 
on this angle. 

We consider a concrete Milne- type problem. Let a dispersive medium 
occupy the half-space x > and let there be a temperature gradient 
along the a:-axis in the medium far from the boundary. We assume 
that radiation at equilibrium with the temperature Tq comes from the 
half-space x < into the medium and the total radiation outgoing 
from the medium is absorbed somewhere and does not come back. 
The temperature far from the boundary has the form 



T = Ti + Kx, 



(13) 



8 



where K is the temperature gradient. Our aim is to find the rela- 
tionship between Tq, Ti, and K. This can be written in the hnear 
approximation as 

Ti - To = fK, (14) 

where / is an unknown coefficient depending on the properties of the 
medium. 

In application to the cosmic gaseous dust cloud, this statement of 
the problem means the following. The radiation at equilibrium with 
the temperature Tq close to 3K is incident on the boundary of a gaseous 
dust cloud. Under cooling of the cloud, energy is transformed and 
radiated outward. Therefore, there is a temperature gradient in the 
cloud. We need to find the relationship between this gradient, which 
is defined by the heat fiow, the temperature Tq, and the temperature 
at the boundary. This relationship is given by formula (14) with an 
unknown coefficient / that must be determined. 

It is easy to show that the following two discrete modes are solutions 
to Eq. (12) (see Sec. 1 below): 

V 

ip+{x,v,uj) = 1, ip.{x,v,uj) = X -. 

The boundary conditions are as follows: 

(f{0,v,u) = 0, 0<v<l, (15) 

(f{x,v,uj) = Ko-\-K[x-—)-\-o{l), X ^ +00, -1 < < 0, (16) 

where K is the temperature gradient (see (13)) and Kq = fKT^^ is an 
unknown coefficient proportional to the relative temperature gradient. 
Note that a similar question for the Boltzmann gas can be stated as 
the Smolukhowski problem on the discontinuity of temperature. 



1. Eigenfunctions corresponding to continuous and discrete 

spectra 
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It can be clearly seen from Eq. (12) and boundary conditions (15), 
(16) that the desired function depends on two variables, x and ii = 
vuj^'^ . Let us perform the change of variables v = jiu'^ and u = u' in 
Eq. (12). The corresponding Jacobian is equal to a;'". Consequently, 
Eq. (12) can be rewritten in the new variables as 

A^ll + ifix, f,) = J uj'^+'E{u)du J ifix, (1-1) 

-l/w" 

Accordingly, under the above transformation, the boundary condi- 
tions (15), (16) for Eq. (1.1) take the form 

(^(0,yu) = 0, 0<//<oo, (1.2) 

ip{x, /i) = Kq -\- K{x — jj.) -\- o{l), a: — )• +00, — oo < yti < 0. (1.3) 
On replacing /i' by u~^ii' in Eq. (1.1), we can rewrite it as 

/^^<^(^,/^) + ^{x,ii) = 

00 1 

= J u^-''E{u)du j ip{x,u-yW- (1-4) 

-1 

Let us perform the separation of variables in Eq. (1.4) as follows: 

X 

(fnix,ii) =e^p{ — )^{r],ii), r] e C. 
As a result, we derive the characteristic equation 

oo 1 

(77 - /i)$(77, /i) = ^^77 / uj''^'E{uj)duj J $(77, u-yW- (1-5) 

-1 

We denote 

00 1 
n(77) = J uj''+^E{uj)duj j ^{ri,uj-"ii')dii' (1.6) 
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and rewrite Eq. (1.5) in the form 

(77 - //)$(ry, 11) = 77T^r/n(77). (1.7) 

2lo{a) 

Equation (1.7) implies that the eigenfunctions corresponding to the 
continuous spectrum are the distributions [11] 

Hv, /^) = 7^7^^^(^)^^^ + 9{v)S{v - /^), ^, e (-00, +00). 

(1.8) 

Here P symbolizes the principal value of the integral, 6{x) is the 
Dirac delta function, and g{r]) can be found from the normalization 
condition for (1.6). Substituting (1.8) into (1.6), we obtain 

n{T])X{T]) = g{T])^ai'n): (1-9) 

where 

1/77° 





and the dispersion function 



CXD 



of the problem can be expressed via the Case dispersion function [12] 

1 



Ac(^) = 1 + / 
-1 

in the following way 



ji — UJ^Z 



00 



A(z) = -i- / J'''^^E{uj)\c{u''z)du. (1.10) 

k[(y-) J 



The boundary values of the function Xc{uj°'z) on the edges of the 

slit 
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are determined by the Sokhotski formulas, 

A±(a;» = J 2 ^' ^ (l.n) 

t Ac(a;», // ^ A,. 

In view of (1.10) and (1.11), the boundary values of X{z) are related 
by the formula 

A±(/i) = A(/i)±z^^, /iGM. 

With the help of expression (1.9), we bring formula (1.8) for the 
eigenfunctions to the form $(77, /i) = $(77, /i)n(?7), where 

Hi, m) = T^rr^^P^— + 7t\^(i - ^) 

2lo{a) T]-/! ^aVn) 
is the eigenfunction corresponding to the normalization condition 

72(77) = 1- 

It can be shown that \{z) has a second-order zero at 2; = 00 
with which two solutions to the original equation (1.4) are associated, 
namely, 

(/?+(a;, //) = 1 and ip-{x, /i) = x — /i. 

2. Eigenfunction expansion for the solution 

Theorem. The boundary value problem (1.2)- (1.4) has a unique 
solution, whose expansion in terms of the eigenfunctions of the char- 
acteristic equation has the form 

00 

(p{x, ii) = Ko + K{x - lA + j e'^'/^Hr], ii)n{r])dr]. (2.1) 



Here the unknowns are the coefficient Kq associated with the dis- 
crete spectrum and the function n{r]), which is the coefficient corre- 
sponding to the continuous spectrum. The expansion (2.1) is inter- 
preted in the classical sense, namely, 

(p{x, 11) = Kq + K{x - 11)+ 
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oo 



1 + M.e-^l.0^(^)^ (2.2) 





where 0+{ii) is the characteristic function of the positive half-Une, i.e., 
e+{ii) = 1 if ;U G (0, +oo) and e+{ii) = if ^ (0, +oo). 

Proof. Using boundary condition (1.2), we pass from expansion 
(2.2) to the following singular integral equation with the Cauchy kernel 
[13]: 

00 

Ko-K., + ^J^-^ + P;^M,) = 0, M>0. (2.3) 



Let us introduce the auxiliary function 

oo 

J r]- z 



Multiplying both sides of Eq. (2.3) by 

^ ■[A+(//)-A-(//)] = 27ri//e,(//) 



(2.4) and using the boundary values N^{/i) and A^(//), we reduce Eq. 
(2.3) to the Riemann boundary problem 

= X-ifi)[N-ifi) + 2loia)iKo - K,fi)], // > 0. (2.5) 
The coefficient in the boundary problem (2.5) is the function 

Consider the following factorization problem for the coefficient G(/i) 



// > 0. (2.6) 
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To solve this problem, we calculate the index x of the coefficient 

argG'(Ai) . 

ZTT L J (0,+oo) 

where the expression [argG(/i)](o,+oo) denotes the increment of the 
function in the brackets under variation of the argument from to 
+00. 

Let 9a[lj) = argA+(//) be the principal value of the argument of 



2 J u''+^E{u)Xc{uj''T)dT 



(^a{lA — arctan — ^ 



00 



Since A+(//) = A (/i) (where the bar symbolizes complex conjuga- 
tion), we have G{ii) = exp(— 2Q,()Li)), and, therefore, 

TT TT 

As a solution to problem (2.6), we take the function 

X{z) = -expy(z), 
which does not vanish at zero, where 

00 

1 f our) - .^^^ 



7T J T 



— Z 


We now return to the problem (2.5) and apply (2.6) to transform it into 
the problem of determining an analytic function from the following 
condition 

X+{ii)[N+{ii) + 2lo{a){Ko - Kni)] = 
+X-(/i)[iV-(/i) + 2loia)iKo - Kill)], fi > 0. (2.7) 
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Taking into account the behavior of the functions entering (2.7), 
we derive the general solution to this problem in the form 

Niz) = -2loia)iKo - Kx) + (2.8) 

where Co is an arbitrary constant. 

In order to use the function (2.8) as the auxiliary function N[z) 
introduced by the formula (2.4), we remove the pole of the solution 
(2.8) at z = oo and make the limit of (2.8) at the point z = oo equal 
to zero by setting 

Co = -2/o(c^)i^, i^o = Vi{a)K, (2.9) 

00 

Vi{a) = -lj iW - M dli. (2.10) 



The coefficient n{r]) for the continuous spectrum can be found from 
Eq. (2.8), namely, 

2«„nW = -2;„(a)if (^-^ - ^ j . 

Thus, the coefficients in (2.1) have been uniquely defined. The 
uniqueness of the expansion (2.1) is proved by the fact that there can 
be no nontrivial expansion of zero with respect to the eigenfunctions of 
the characteristic equation. Obviously, expansion (2.1) automatically 
satisfies to the boundary condition (1.3) and (1.2) holds by construc- 
tion. Direct verification shows that the expansion (2.1) satisfies to the 
Eq. (1.4). The theorem is proved. 

3. Numerical results. Discussion and conclusion 

First, we note that we can find the exact value of the integral (2.10) 
for a = 0. Indeed, if a = 0, then it can be seen from the formula (1.10) 
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that the dispersion function X{z) coincides with the Case dispersion 
function Xc{z). In this case, 

oo 

1 f TTT 

= — / arccot— —-c/r = 0.71045. (3.1) 
TT J 2Ac(rj 



Let us find an approximate value of the integral (2.10) for a > 0. 
Applying the saddle-point method, we replace X{z) by an approximat- 
ing function X{z). The expression (1.10) implies that the saddle point 
ujQ is determined by the equation 

^ ^ ^ (3.2) 

The approximate solution Cjq to Eq. (3.2) for « + 4 > 1 can be 
calculated by the formula 

a;o = (a-t-4)(l-2e-"-4). (3.3) 

Note that formula (3.3) provides a satisfactory approximation even 
for a = 0. Indeed, in this case ujq = 3.83002 and a;o=3. 85347. For 
a = 2, we have a;o=5. 96941 and cDq = 5.97025, i.e., formula (3.3) gives 
a nearly exact value of the saddle point. According to the saddle-point 
method, we have 

A(^) = Xciuj^z) 

and, therefore. 



Vi = -- I [eair] - 7r]dT = -- I arccot-^^^ 





— J l^aiT) — ^]dT = — J arccot - 



whence 

Vi=uj,^V,\ (3.4) 

Note that formula (3.4) with A = exactly implies the formula 
(3.1) and for a = 2, we have Vi = 0.01994. The decrease in Vi 
with increasing a relates to the fact that the contribution from high 
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frequencies (large values of cu) is reduced during the transport process 
in the dispersive medium due to the faster growth of the collision 
frequency iy{co)) with increasing a. At the same time, the radiation 
from the boundary surface is not limited by this factor nor does the 
high-frequency contribution depend on a. 

The result is that the relative outflow of energy from the surface 
increases together with a (as compared to the volumetric outflow). 
As a result, the coefficient Vi, which characterizes the temperature 
drop between the boundary of the exterior medium and the incident 
radiation, decreases. 

In conclusion, we note that in the present paper, a generalization 
of the Case — Zweifel classical theory [12] to the transport equation 
for a massless Bose gas was considered for the first time. It turns out 
that the suggested method permits transport equations with a double 
integral on the right-hand side to be solved analytically. 

This work is supported by the Russian Foundation for Basic Re- 
search, Grant No. 97-01-00333. 
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